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ABSTRACT 

We report 850 /mi Sub-millimeter Array (SMA) observations of four gravitationally lensed submil- 
limeter galaxies (SMGs), A370-2, A2390-1, A2390-3 and A2390-4, which were originally discovered 
with the Submillimeter Common-User Bolometer Array (SCUBA). Our SMA detection of A370-2 with 
a submillimeter flux of 7.95 ± 0.60 mJy unambiguously identifies the counterparts to this source at 
optical and radio wavelengths. A2390-1 is an ultraluminous infrared galaxies with a submillimeter flux 
of 5.55 ± 0.92 mJy and a redshift of 1.8 ± 0.2 computed from submillimeter/radio flux ratio analysis. 
We resolve A2390-3 into two components, A2390-3a and A2390-3b, with fluxes of 3.15 ± 0.63 mJy 
and 1.92 ± 0.60 mJy, respectively. The structure of the system could be consistent with morphologi- 
cal distortion by gravitational lensing. The lack of counterparts in the optical and infrared indicates 
a heavily dust-enshrouded nature, and a non-detection in the radio implies that these two sources 
probably lie at z > 4.7, which would make them among the most distant SMGs known to date. Our 
non-detection of A2390-4 suggests either that there are multiple fainter submillimeter sources within 
the SCUBA beam or that the SCUBA detection may have been false. Our precise positions allow 
us to determine accurate amplifications and fluxes for all our detected sources. Our new results give 
a shallower power-law fit (-1.10) to the faint-end 850 /im cumulative number counts than previous 
work. We emphasize the need for high-resolution observations of single-dish detected SMGs in order 
to measure accurately the faint end of the 850 /im counts. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation — galaxies: 
high-redshift — submillimeter: galaxies 



1. INTRODUCTION 

During the past decade, deep blank-field submillime- 
ter/ millimeter surveys from the Submillimeter Common- 
User Bolometer Array (SCUBA) on the James Clerk 
Maxwell TelescopeF] (JCMT) and the Max Planck Mil- 
limeter Bolometer (MAMBO) array on the IRAM 30 me- 
ter telescope have resolved sources brighter than ~ 2 mJy 
that account for ~ 20%- 30% of the 850 /xm extragalactic 
background light (EBL) (|Barger et al.|1999a[ Eales et al. 
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The rest of the 850 /tm EBL can be accounted 
tor by sources over t he 0.3-2 mJy range measured using 



gravitational lensing (Blain et al. 1999 Cowie et al. 2002 
Knudsen etldj]2008| ). 

However, further studies on individual submillimeter 
galaxies are hampered by the relatively poor resolution 
of single-dish submillimeter telescopes (e.g., 15" FWHM 
at 850 /im on SCUBA). This causes a correspondingly 
large uncertainty in the source positions and confusion 
in identifying counterparts to the submillimeter sources 
at other wavelengths. It also makes the gravitational 
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lensing uncertain, particularly for sources with very high 
amplification. Techniques have been developed to work 
around this problem: The empirical correlation between 
non-thermal radio emission and thermal dust emission 



(Condon 1992L together with the high astrometric pre- 



cision ot radio interfero meters, has been widely used as 
a tool to study SMGs. Chapman et al. (20051 deter- 
mined the redshifts of 73 radio-identified SlVlGs using 
Keck spectroscopy, finding a median redshift of ~ 2.2. 
It has also been found that SMGs appear to have simi- 
lar properties to the local ultraluminous infrared galax- 
ies with far-infrared (FIR) lu minosities of ~ 10 12 -10 13 
L Q ( |Sanders fc Mirabel||l99"6"l ). The lack of strong X-ray 
emission from submillimeter sources (Lx/Lfir = 0.004) 
suggests that most of the FIR/submillimeter flux is dom- 
inated by dust-reradiated emissio n from intense star for - 
mation rather than AGN output ( Alexander et al.| 2005), 
which makes studying the FIR/submillimeter sources a 
path to robustly determining the dust-obscured star for- 
mation history. 

There are several drawbacks to using radio emission to 
identify the SMGs. Firstly, with the current sen sitivity (5 



a ~ 20 /xJy at 20 cm) of radio interferometers ( | Morrison 
et al.|20iu ), the radio-identified SMGs are mostly bright 
(» 2 mJy at 850 /im) , so their properties may not be rep- 
resentative of the submillimeter population as a whole. 
Secondly, the radio flux drops at high redshift due to 
the positive if-correction, whereas the submillimeter re- 
mains almost invar iant over the redshift range z ~ 1 — 8 
(Blain et al. 2002). Thus, the radio-identification tech- 



nique is biased against high-redshift sources. Lastly, 
there are cases where multiple submillimeter sources are 
located within the beam of a single-dish submillimeter 
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telescope (Wang et al. 20111, so a radio source within 
the beam may not be the correct or only counterpart. 

In the last few years our ability to address these issues 
has greatly improved with the advent of te lescop es like 
the Submillimeter ArrajQ (SMA; |Ho et al.|[2004| . The 
SMA provides imagi ng in the sub millimeter regim e with 



arcsecond resolution ( Iono et al.| :2006 W ang et al.||2007[ 
Younger et al.||2007| |2UolT 
et al. 2010 ) 
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Some distances to submillimeter sources 



have been found from spectro scopic redshifts (Capak 
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strong radio coun terparts \ Younger et al. 2007 2008 
Cowie et all 2009). This again illustrates the problems 
with the radio identification technique. However, most 
of the sources studied so far have fluxes around 10 mJy, 
and they constitute only a small fraction of the submil- 
limeter EBL (~ 2%, according to the counts in Coppin 
et al 2006). The nature of more typical SMGs with sub- 
millimeter fluxes < 2 mJy is still an open question. 

With the improved sensitivity from the recently up- 
graded double bandwidth (4 GHz) on the SMA, we are 
now able to reach a sub-mJy level of sensitivity and probe 
relatively submillimctcr-dim sources. In order to sample 
very faint sources we have used this capability to ob- 
serve sources in the gravitationally l ensed regions of tw o 
massive clusters, A370 and A2390 ( |Cowie et aLl|2002[ ), 
which enables us to take advantage of the lensing am- 
plifications. We list the basic information on our target 
sources m Table [TJ where we follow the index sequence 
from the catalog of |Cowie et al (2002). We give the 
observational details and discuss our reductions in Sec- 
tion 2. In Section 3 we report on our SMA observations 
of one radio-bright source, A370-2, and two radio-dim 
sources, A2390-3 and A2390-4. We show unambiguously 
whether there are counterparts to these sources in the 
optical, infrared, or radio. In Section 4 we discuss the 
implications of our new measurements on the faint-end 
850 /im galaxy number counts. We summarize our results 
in Section 5. We assume the WMAP (Wilkinson Mi- 
crowave Anisotropy Probe) cosmology throug hout: Hp 
= 70.5 km s" 1 Mpc 
et~aLl[2009 ). 
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2. OBSERVATIONS AND DATA REDUCTION 

Seven 345 GHz SMA tracks of A370-2, four tracks of 
A2390-1, three tracks of A2390-3 and three tracks of 
A2390-4 were taken between 2007 and 2010. We used the 
subcompact configuration (16.5 - 32 meter baselines) on 
two tracks of A370-2 and the compact configuration (20 - 
70 meter baselines) on all other tracks in order to obtain 
a better signal to noise ratio. The shape of the band- 
pass was measured by observing two quasars 3C454.3 
and 3C111. The sources 0238+166, 0132-169, 2148+069 
and 2232+117 were observed for time dependent gain 
calibration, and the planets Uranus, Ceres and Callisto 
were used as flux calibrators. We used the data reduc- 
tion package MIR to calibrate the visibilities, and we 

6 The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics and is funded by the Smith- 
sonian Institution and the Academia Sinica. 



produced t he images using the MIRIAD routines (Sault 
eFal][l995 



To minimize the noise level we combined all the tracks 
for each source and made the final continuum images 
with naturally weighted baselines. The cleaned images 
(the outcome of the C1EAN algorithm) of the calibra- 
tors from each and every track show only one cleaned 
component at the phase center of the calibrators. After 
performing RESTOR to convolve the cleaned component 
with the point spread function, every calibrator appeared 
as a point source, which is an indication of good phase 
calibration. In order to prevent over-cleaning of the dirty 
maps, we also made sure that the r.m.s of the cleaned 
maps, excluding the detected sources, agreed with the 
theoretical noise. One should keep in mind that there is 
a typical ~ 10% error introduced by the flux calibration, 
on top of the image r.m.s noise for flux measurements. 

The theoretical r.m.s noise at the phase center and the 
resolution of each of the final images of the four SMGs 
are given in Table [T] Note that the better weather con- 
ditions for A2390-3 and the double bandwidth (4 GHz) 
used in the last track of the A2390-3 data make the sen- 
sitivity comparable to A370-2, even though there were 
four fewer tracks of data for A2390-3 than for A370-2. 
The sensitivity of the image of A2390-1 is not as good as 
others due to bad weather conditions during the first two 
tracks of data. For A2390-4 all three tracks of data are 
double bandwidth, so a better sensitivity was obtained. 

3. RESULTS 

We have a strong detection (> 5 a) toward A370-2 
and A2390-1, respectively, and two significant detections 
(both > 3 a) toward A2390-3. However, we found no 
sources at the A2390-4 position, even with these deepest 
SMA observations toward a single SMG. It is possible 
that there is more than one source contributing to the 
submillimeter emission in this region, and that once we 
are able to resolve it into individual sources, their fluxes 
will be found to be below the current detection limit, 
as has been the c ase for several su bmillimeter sources 
in the GOODS-N ( |Wang et al.|20lT] >. Alternatively, the 
SCUBA detection may have been false. Again, this is one 
of the main reasons to use high-resolution telescopes to 
observe SCUBA sources. We discuss the A370-2, A2390- 

I and A2390-3 detections below. 

3.1. A370-2 

A370-2, which is a l so referred to as S MM J02399-0134 
(|BargCTet^l.|1999bllSm"ail et al.|[2000| or SMM J02396- 
0134 ( |Greve et al ||2005| |Wu et al| |2009), has three op- 
tical sources within the SCUBA beam (Figure [lj). The 
galaxy 13 has a distorted ring morphology. Its formation 
mechanism has been suggested to be either the dynami- 
cal axial penetration of a smaller companion galaxy into 
a disk galaxy or a pair of superbubb les driven by the 
intense starbu rst at the central region ( Taniguchi & Mu- 
|rayama||2001[ ). A spectroscopic redshift of ~1.06 was in- 
dependently obtained by two groups using IRIS on Keck 

II and OSIS-V on the Canada-France-Hawaii Telescope 
(CFHT), respectively, and the identification of the high- 
ionization lines of [ Ne V] indicates that 13 is likely to 



be a Seyfert galaxy (Barger et al. 1999b Soucail et al 



19991. located beyond the cluster A370, the signal from 



13 is amplified by a factor of 2.3 due to gravitational 
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TABLE 1 

SMA OBSERVATIONS 



Target 


R.A. (2000) 


Dccl. (2000) 


SCUBA Sg 50 ^ m 


SMA S850 |un 


Amplification 


Theoretical Noise (1 a) b 


Resolution 6 


Source 


(h m s) 


(d m s) 


(mjy) 


(mjy) 


(min,max) 


(mjy beam -1 ) 




A370-2 


02 39 56.63 


-1 34 27.0 


6.68 ± 0.58 


7.95 ± 0.60 


2.3 


0.57 


2'.'4 x 2'.'3 


A2390-1 


21 53 33.58 


17 42 42.3 


7.57 ± 0.93 


5.55 ± 0.92 


1.9 (1.8,2.0) 


0.77 


2'.' 3 x 179 


A2390-3 


21 53 35.48 


17 41 09.3 


3.24 ± 0.78 


5.07 ± 0.87 


52 (0.6,52) 


0.55 


2'/0 x 1'.'6 


A2390-4 


21 53 38.21 


17 41 52.3 


2.64 ± 0.72 


< 0.96 c 


11 (> 6.7) 


0.32 


2'.' 2 x 1'.'8 



" IjNS'JXXJ SCUBA 850 ^m fluxes and source amplifications, which were determined using the SCUBA positions and 

b Obtained from natural weighted baselines. The resolution is the FWHM of the synthesized beam. 
c The 3 a upper limit for a point source. 




Fig. 1. — Multiwavelength images of the ring galaxy A370-2 and its surroundings: Left: HST ACS F625W image - Center: SMA 850 fim 
contours overplotted on the F625W image - Right: 20 cm VLA contours overplotted on the F625W image. Note that all three images are 
centered at the position given in Table [l] The big black circle in each image shows the SCUBA beam size (15" X 15"). North is up and 
East is to the left. 



lensing ( Cowie et al. 2002 ) . A stron g CO molecu l ar de- 



2005|) 

4Tand 



tection toward L3 was reported by |Greve et al 
using the IRAM Interferometer. L4 is at z = 
appears to be a background normal galaxy. The bright 
passive elliptical gal axy L5 is one of the cluster members 
at z = 0.37 (#32 in |Mellier et"aL|l988 ). 

The evidence discussed above implies that L3 is the 
only active source in this region, and it has been argued 
that L3 is t he likely counterpart of the SMG JBarger 



et al. 



1999b). Our SMA data confirm this identification 



and clearly demonstrate that the submillimeter emission 
is from the center of galaxy L3 (Figure 1). A370-2 is 
detected both in our dirty and CLEANed images. The 
point-source fit routine IMFIT in MIRIAD with a 10" 
box centered at the peak position produces a flux of 
7.95 mjy and also yields an r.m.s for the residual image 
of 0.60 mJy/beam. Note the flux measurements through- 
out this paper are all primary beam corrected. The fitted 
peak position is at a(J2000.0) = 2' l 39 m 56.55 s , <5(J2000.0) 
= -1°34'26.54" with an error of 0715. 
The total submm flux of 7.95 ± 0.60 mj y agrees with 



the SC UBA flux of 6.68 ± 0.58 mjy from |Cowie et al. 
( 2002 ) (Table [I]) . We u se the spectral energ y distribu- 



tion (SED) model from |Barger et aL (2000) based on 
an Arp 220 template with an assumed dust temperature 
and extinction coefficient (Tj = 47 K, /3 = 1) to com- 
pute the infrared luminosity. The SED in the infrared 



and radio regime is scaled to fit our 850 /im flux, and the 
infrared luminosity can simply be calculated from the 
area under the SED curve between 8 /im and 1000 /im. 
With the adopted redshift of 1.06, the de-lensed flux 
of 3.46 mjy corresponds to an infrared luminosity of 
~ 2.5 x 10 12 Lq, which is typical of ultraluminous in- 
frared galaxies (ULIRG, L IR > 10 12 L Q ). 

A strong 1.4 GHz (20 cm) detection of the source was 
repo rted with ~ 570 r esolution using the Very Large Ar- 
ray ( Smail et al.|2000 ). The emission peaks at a location 
close to L5, and the morphology appears unresolved and 
elongated from north-east to south-west. In order to re- 
solve this region, we made use of the higher resolution 
(1768 x 1749) VLA archival data taken on 1999 August 
at the same frequency toward this region in the most ex- 
tended A configuration. The data clearly show two radio 
sources at the locations of the optical sources L3 and L5 
(Figure 1). Note that without our SMA observations, 
it would still be ambiguous which radio source is the 
counterpart of A370-2 based on the correlation between 
submillimeter and radio emission. 

The radio and submillimeter peak positions at the L3 
location are still offset by 076, which may be caused by 
the contamination of the jet structure extending from 
L5. The radio flux from L3, 0.573 mjy, is estimated 
from a 275 radius circular aperture covering most of the 
radio flux from L3. We ran a Monte Carlo simulation 
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Fig. 2. — A2390-3 SMA image showing the two submillimetcr 
detections marked by the letters a and b. The black contours 
have levels of —2, 2, 3, 4, 5 a, where 1 a is 0.58 mjy beam -1 . The 
gray scale image shows the linear scale range from 1.0 to 5.0 mjy 
beam -1 . A 1"9 X 1'.'6 beam is presented in the lower left corner. 
The green ellipse is the outcome of gravitational lensing, assuming 
a 1" radius source in the middle of two sources on the source plane. 

with the same aperture on the VLA image, and this 
gave a noise level of 0.011 mjy. Together with the sub- 
millimeter flux, our results show a S^o A , m to Si. 4 qhz 
ratio of 12.6 - 15.2. For sources lying at z < 3 in 



the R ayleigh- Jeans long- wavelength limit, |Barger et al. 
(2000) obtained a formula describing the relation be- 
tween the redshift and the submillimcter/radio flux ratio, 
z = 0.98(5850 ^m/Si.4 ghz)°' 26 -1- This implies z = 0.89 
- 0.99 for L3, which matches well to the measured spec- 
troscopic redshift of z — 1.06. Generally, AGN contam- 
ination could be one o f the reasons responsible for the 
slight underestimation ( Carilli fc Yun|[l999 Wang et al. 
20071. In this particular case, however, an inaccurate 



assumed dust temperature and radio flux contamination 
from the extended jet structure could also bias the esti- 
mation downward. 

3.2. A2390-1 



Cowie et al | (2002) discovered the bright SMG A2390- 
m their SCUBA survey, and its 850 /zm flux is 7.57 ± 



0.93 mjy. |Metcalfe e t al. (2003) performed a deep survey 
on A2390 using the 1SOCAM on ESA's Infrared Space 
Observatory (ISO) and found a ~ 3.7 a detection at 15 
fim toward A2390-1. Our SMA data confirm the SCUBA 
detection with a flux of 5.55 ± 0.92 mjy obtained from 
the point-source fitting routine IMFIT with a box size 
of 10". The fitted position from IMFIT is a(J2000.0) = 
2l' i 53 m 33.33 s , c5(J2000.0) = 17°42'50.05" with an error 
of Of! 7. Note that we only use the last two tracks of data 
to determine the source position because of high phase 
noise in the first two tracks of data. 

There is a point source-like radio counterpart of A2390- 
1 in our latest A2390 VLA image (detailed information in 
Section 3.3). We compute the peak position and the flux 
by fitting a two dimensional gaussian with fixed major 
axis, minor axis, and phase angle obtained from the VLA 



synthesized beam. The radio position is at a(J2000.0) 
= 2l' l 53 m 33.31 s , <5(J2000.0) = 17°42'50.33" with 0'.'03 
error, which agrees with our SMA position. We adopt 
the position of A2390-1 measured from the VLA image, 
since it has a better accuracy. We use the VLA position 
and LENSTOOL to determine an amplification of 1.8 
for A2390-1, which agrees well w ith the SCUBA mea- 
surement from |Cowie et al" (2002). 

The fitting algorithm also gives a radio flux of 111.5 ± 
6.5 /zJy (Table [2]), which allows us to obtain the redshift 
of A2390-1, from the submillimeter/radio flux ratio tech- 
nique. We find z = 1.8 ± 0.2. The infrared luminosity of 
A2390-1 that we estimate from the SED model described 
in Section 3.1 is 2.9 x 10 12 L Q . This shows that A2390-1 
is also a ULIRG. 

3.3. A2390-3 



A2390-3 was first identified by |Cowie et al.| ( [2002 1 
in their SCUBA map with a total flux of 3.24 ± 
0.78 mjy. With seven times higher resolution, our SMA 
map detects the submillimeter object and appears to 
resolve it into two sources (Figure [2]), which we re- 
fer to as A2390-3a and A2390-3b. The simultaneous 
point-source fitting routine IMFIT with a box size of 
10" gives fluxes of 3.15 mjy and 1.92 mjy, respectively. 
The r.m.s of the residual image is 0.58 mjy beam -1 , 
whereas the primary beam corrected noise is 0.63 mjy 
for A2390-3a and 0.60 mjy for A2390-3b. The fitted 
positions are a(J2000.0) = 2l' i 53 m 35.141 s , <5(J2000.0) = 
17°41'05.9" with a 074 error for A2390-3a and a(J2000.0) 
= 21 /l 53 m 35.21 s , <5(J2000.0) = 17°41'07.5" with a 076 
error for A2390-3b (Table m. The total flux of 5.07 ± 
0.87 mjy is consistent witn the measured SCUBA flux. 
The error of 0.87 mjy is obtained from error propagation. 

We optimized the output image (Figur e [2]) by apply- 
ing the robust weighting of Briggs ( 1995 ) with a robust 
parameter of +1.0 for the sake of making a balance be- 
tween resolution and signal-to-noise ratio (S/N). To test 
the robustness of the detections, we iterated the pro- 
cess of making the image with different weightings on 
the Fourier transformation and different box sizes for 
the CLEAN algorithm. Noise spikes would vary or even 
disappear with the iterations; however, the signals from 
A2390-3a and A2390-3b remain, and the detections are 
always significant (> 3 a for each compon ent). 

The amplification factor obtained by |Cowie et al.| 
( 2002 ) on this source varies from 0.6 to 52 due to the posi- 
tional uncertainty and redshift indeterminacy. However, 
with our arcsecond level precision, we can now obtain 
more accurate amplifications using LENSTOOL, which 
models the effects of the gravitational lensing by taking 
three dimensional mass distributions within the cluster 



into account ( Kneib et al. 1996 1 



We took the fitted positions and assumed that A2390- 
3a and A2390-3b are both point sources with a 1" ra- 
dius size on the image plane. LENSTOOL calculates the 
positions, semi-major axes, semi-minor axes, and incli- 
nation angles of the sources on the source plane. The 
amplification is then obtained by taking the inverse of 
the product of the semi-major and semi- minor axes (Ta- 
ble [2]) . We assume the source plane is located at z — 5 
for the LENSTOOL calculation. The amplifications of 
sources lying beyond z = 1 with modest amplifications 



ACS F850LP 




Fig. 3. — Multiwavelength images of A2390-3. North is up and East is to the left. The size of each image is 30" X 30". The big white 
circle in each image is centered at the position given in Tableland represents the SCUBA beam size (15" X 15"). The SMA position of 
A2390-3a is labeled in each image with a 2''0 diameter yellow circle. The gray scale images have inverse scales. Three of the IRAC channels 
(1, 2, and 4), corresponding to 3.6 (im, 4.5 (im, and 8.0 fim, are presented in the combined IRAC image with the color codes labeled. 



have little dependence on redshift |Blain et al. (1999 1 . 
We find an amplification of 2.3 for A2390-3a and 2.9 tor 
A2390-3b. 

We also present archival images from the Hubble Space 
Telescope (HST) Advanced Camera for Surveys (ACS) 
z-band (F850LP), the Spitzer Infrared Array Cam- 
era (IRAC) 3.6 /im, 4.5 /im, and 8.0 /xm bands, and 
the Multiband Imaging Photometer for Spitzer (MIPS) 
24 /an band (Figure p| . We present radio maps from 
VLA data obtained on October 20, 2008 at 1.4 GHz 
(Wang et al., in preparation), as well. With the most 
extended A configuration, the 1.4 GHz synthesized beam 
has a size of l'/4 x l'/4, and the 1 a noise level is 6.5 /zJy 
beam -1 . Two cluster members lie within the SCUBA 



beam: CNOC 100890 was discovered by |Yee et ail ( 1996) 
from the CNOC cluster redshift survey using the CVriT, 
while 0657.7+0705.0 was reported from a strong Ha 
emission line search ( |Balogh fc Morris 20001. A strong 
amplification toward this region could be caused by these 
two galaxies. 

Interestingly, Figure [3] shows that there is no apparent 
detection corresponding to the SMA position (the SMA 
position of A2390-3a is shown) at other wavelengths. 
This evidence implies that both sources must be heav- 
ily embedded in dusty envelopes and/or at very high 
redshifts. For each source, the 1.4 GHz flux was mea- 
sured from a given circular aperture which covers most 
of its submillimeter flux (Table [2| , yet it is still a non- 
detection. We assigned zero to the VLA flux and adopted 
the 1 o~ errors obtained from Monte Carlo simulations as 
the upper limits on the radio fluxes, which are 5.2 fiJy 
and 3.3 /iJy for A2390-3a and A2390-3b, respectively 

Note that we are aware of the possibility of contamina- 
tion from the residual sidelobes of the strong radio-bright 



cD galaxy in A2390. However, with the help of the latest 
data reduction techniques provided in AIPS, we find that 
the residual sidelobes from the cD galaxy can be mostly 
removed, except for the area close to the cD galaxy and 
lying generally N-S and away from the source. A2390- 
3a/b are located far from that area. Also, in our Monte 
Carlo simulations we avoid regions affected by bright 
sources and conspicuous residual sidelobes. Thus, we 
believe that our estimate of the 1 a flux limit is robust. 

With the absence of a detection at all wavelengths ex- 
cept the submillimeter, it is not possible to obtain op- 
tical or infrared spectroscopic redshifts or photometric 
redshifts. A natural option would be to use the submil- 
limeter/radio flux ratio analysis. Since the equation in 
Section 3.1 is limited to lower redshift sources (z < 3), 
.eral equations (Eqs. (2) and (4) in 
. We find z > 5.2 for A2390-3a and 
3b. If confirmed with an identifica- 
tion of CO lines, they could be among the most distant 
SMGs known. However, the high submillimcter/radio 
flux ratios could also be caused by using an incorrect dust 
temperature. In this case, the non-detection in the mid- 
infrared is a tentative indication that the dust tempera- 
ture is likely to be low. Also, SMGs tend to have lower 
dust temperatures than their local infrared counterpart s 
( [Pope et al.|2006)|Huynh et al |2007[ [Hwang et al.|2010 ) 



Barger et al. 


2000 


z > 4.7 for 1 


12390 



Any dust temperature lower than our adopted one would 
move the source even further (i.e., to higher redshifts). 
Thus, our estimation of the redshifts should be reason- 
able. 

Having the lower limit redshift and amplification in- 
formation for both sources in hand, we can also com- 
pute lower limits to the infrared luminosities (Table [2]). 
The star formation rates can be estimated if we as- 



G 



TABLE 2 

Detected sources in the SMA observations 



Source 


R.A. (2000) 


Decl. (2000) 


S850^m 


&1.4GH* 


z" 


R.A. Offset 


Dec. Offset 


Amplification 


Lir 




(h m s) 


(d m s) 


(mjy) 


(MJy) 




(arcsec) 


(arcsec) 






A370-2 


02 39 56.55 


-1 34 26.5 


7.95 ± 0.60 


573 ± 11 


0.89 - 0.99 


-1.2 


0.5 


2.3 


2.5 x 10 12 


A2390-1 


21 53 33.31 


17 42 50.3 


5.55 ± 0.93 


111 ± 7 


1.60 - 1.99 


0.0 C 


1.0 C 


1.8 


2.9 x 10 12 


A2390-3 






5.07 ± 0.87 


-4.0 ± 6.2 


> 6.2 






2.5 6 


> 1.7 x 10 12 


A2390-3a 


21 53 35.14 


17 41 05.9 


3.15 ± 0.63 


-4.5 ± 5.2 


> 5.2 


-4.8 


-3.4 


2.3 


> 1.1 x 10 12 


A2390-3b 


21 53 35.21 


17 41 07.5 


1.92 ± 0.60 


0.5 ± 3.3 


> 4.7 


-3.8 


-1.9 


2.9 


> 4.6 x 10 u 



b Estimated under the assumption of a point source located ^ 
c The SMA pointin g position is different than the one reported in 



Cowie et 



lAp. 



2002 



middle of A2390-3 


Cowie et al. 


2002; 



Here we use the reference position from 



sume starburst galaxies and a standard Salpeter ini- 
tial mass function. In thi s work, we use the formula 
M = UxlCT 10 L FIR /L Q ( |Kennicutt|l998| , which gives 
188 and 78 M /yr for A2390-3a and A2390-3b. 

It is possible that these two sources are coincidcntly 
passing through the line-of-sight at different distances. 
Our estimates of the redshifts are rough because of the 
uncertainty in the radio-FIR correlation. LENSTOOL 
shows that the separation between the two sources goes 
from l'/8 on the image plane down to 1"7 on the source 
plane. More interestingly, if we assume a point source 
located in the middle of the two sources on the source 
plane with a 1" radius morphology, then the gravitation- 
ally lensed image on the image plane appears elongated 
and covers both sources (green ellipse in Fig. [2]). This 
implies that it is also possible (perhaps more likely) that 
these two sources are the outcome of the distortion by 
gravitational lensing of a close pair or of one source bro- 
ken into two. At this point we cannot be secure about the 
nature of this object, and more rigorous investigations, 
such as CO line searches, would be needed to proceed 
further. We have also tried treating the object as one 
source, and we present that information in Table [2j as 
well. 

4. DISCUSSION 

Previ ous SMA studies have focused on very luminous 
SMGs (llono et al.|2006"l [Younger et al.|2007| [Wang et al. 
|2007[|Cowie et al.|2009p Those bright SMGs have a total 



infrared luminosity in the range of 10 



12 



10 



13 



de- 



pending on the redshift. While A370-2 and A2390-1 are 
also ULIRGs, compared to those giant galaxies, A2390- 
3b is much more typical, thanks to the gravitational lens- 
ing. Even if we assume that A2390-3b has a redshift of 
6, the infrared luminosity still remains below 5 x 10 11 
Lq. Moreover, our non-detection toward A2390-4 may 
also imply multiple faint sources. Given our detection 
limit (0.96 mjy for 3 a), the assumed redshift (z = 3 
obtained from the submillimeter/radio flux ratio), and 
the amplification (6.7), the total FIR luminosity would 
easily be of the order of typical galaxies (~ 1.4 x 10 11 
Lq). 

Giant SMGs tend to dominate the univers al star for- 
mation hist ory during the ep och z = 1 - 3 ( |Chapman 
et al. ||2005| [Wang et alT||2006[ ). The answer to the ques- 
tion ot whether this situation continues to higher red- 
shifts is critical for understanding the star formation his- 
tory. Smaller galaxies like A2390-3b are expected to 
dominate the light in the very early Universe (z ~ 5 



o 



10 2 
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Fig. 4. — Updated cum ulative 850 nm sou rce counts (filled cir- 
cles) using the method of |Cowie et al.| | |2002| (open squares). The 
solid curve is an area-weighted maximum likelihood power-law fit 
to the updated data points over the range 0.1 to 4.0 mjy with a 
slope of —1.10. The dashed curve is the fit to the original data and 
has a slope of —1.14. The typical error bar shown in the bottom 
left corner assumes Poisson statistics. 



or higher), or else bottom-up clumping cold dark mat- 
ter models of galaxy formation will be seriously chal- 
lenged. Thus, obtaining the redshifts for A2390-3b and 
other submillimeter-faint sources, possibly through CO 
observations, is an important next step. 

With our updated positional information, we have re- 
co mputed the 850 cumulative number c ounts shown 
in|Cowie et al. (2002). We used the data in Cowie et al. 
( |2002 1 and updated the three sources in this work with 
newly obtained information. The results are sh own in 
Figure |i} W e adopt the method described in |Cowie 



et al. ( |2002[ ) to calculate the source areas and cumu- 
lative counts. We calculate the upper and lower error 
bars usi n g Pois son statistics, and we use the table from 
Gehrels (11986) to find the confidence limits on small 
numbers. We also perform an area- weighted maximum 
likelihood fit to the updated data over the range from 0.1 
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to 4.0 mJy, which gives the best fit equation of 
N(> S) = 8.2 x lO 3 ^- 110 . 

The slope of the fit is slightly shallower than the previous 
fit (—1.14) using the same clusters. Note that we assume 
no source in A2390-4 and take A2390-3 to be a single 
source. 
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Fig. 5. — Cumulative 850 fim source counts from other lensing 
surveys compared with the updated source c ounts from this work . 
For clarity we only show the error bars on the Knudsen et al.| | |2008[ l 
data points. The typical 1 a error bar for our counts is shown in 
the bottom left corner. The solid curve is our power-law fit to the 
updated source counts wit h a slope of -1.10 , the da shed curve is the 
double power-law fit from [Knudsen et al. (2008 ). an d the dotted 
curve is the parametric fit trom [riarger et al.| ( |l999a[ | . 



In Figure [5] we compare our updated results at the 
faint end (0. 1 mJy < 5*850 < 4 mJy) with previous lens- 
ing survey s (Bl ain et al. 1 19991 Cowie et al. 120021 iSmail 
et al |2002| |Knudsen et aL|2008| ) . We overplot various fit- 
ted curves: tfie solid line is the power-law fit in this work 
th e dashed line shows the double power-law fit described 



Knudsen et al. (2008), and the dotted curve shows 



an em pirical tit to a blank field survey ( Barger et al 



1999a) that was constrained by the 850 /Ltm EBL mea- 



surements at lower fluxes. All three curves reasonably 
describe the cumulative counts, though Barger et al.'s 
model may be preferred as a better match if the trend of 
flattening continues to counts lower than 0.1 mJy. Our 
faint-end counts are slightly lower but still in a good 
agreement with the previous literature. 

The changes in our derived counts emphasize the un- 
certainty introduced in the lensing analysis when the 
source positions are poorly determined. The reason for 
the change is that the faintest sources get a large boost 
from gravitational lensing (amplifications > 10), and for 
these sources the variations in the amplifications due to 
positional errors can be an orde r of magnitude (Cowie 
et al. 2002 Knudsen et al. 2008). The small source ar- 



eas at the faint fluxes then cause large changes. Using 
A2390-3 as an example, the accurate position from our 
results brings the amplification from 52 down to 2.5 as- 
suming one source. In fact, our updated flux and ampli- 
fication for A2390-3 is the main reason for the shallower 
fit. 
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Fig. 6. — Same as Figure[4] but with horizontal error bar s shown 
on the 850 /im fluxes of the sources in |Cowie et al. (2002 ), which 
are caused by amplification uncertainties. 



In Figure|6]we show the Cowie et al. ( 2002 ) data points 
with error bars on the 850 /zm fluxes, which are due 
to the uncertainties in the amplifications (i.e., resulting 
from the positional uncertainties from the SCUBA mea- 
surements). Large flux errors on the faint-end sources 
clearly demonstrate the problems of using measurements 
from single-dish telescopes. The fact that our shallower 
power-law fit results from our accurate determination of 
the amplification of A2390-3 implies that accurate posi- 
tions of SMGs are not only critical for finding the correct 
counterparts in order to determine the nature of individ- 
ual SMGs, but also for measuring the overall properties 
of the submillimeter sources. 
Moreover, the multiplicity of SMGs has b een implied 



from several in direct studies using 24 /im ( Pope et al 
2006I) 350 urn dKovacs et al.||2010| ), and radio ( jlvisori 
eraT||2007| [Younger et al .||2009C, as well as CO interfer 
ometry (ITacconi et al.||2 



|2O06[ |2008l lEngel et al||2010 > 



The most recent discovery by Wang et al. ( |201ip 



the GOODS-N region that some SCUBA-detectecFSMGs 
split into multiple fainter sources in deep SMA imaging 
supports this result. This multiplicity of SMGs can po- 
tentially increase the number counts at the low luminos- 
ity end and may also fundamentally change the luminos- 
ity function of SMGs. Thus, future high-resolution ob- 
servations using interferometers like the SMA or ALMA 
will be critical for obtaining a better understanding of 
the star formation history. 

5. SUMMARY 



8 



We have reported the results of our SMA observa- 
tions of four SCUBA-detected sources, A370-2, A2390-1, 
A2390-3 and A2390-4. Although it had been suggested 
with indirect evidence that the optical source L3 was 
the likely counterpart of A370-2, our direct detection 
of submillimeter emission from L3 unambiguously con- 
firms this. A2390-1 is a ULIRG with a submillimeter 
flux of 5.55 ± 0.92 mJy and a redshift of 1.8 ± 0.2 com- 
puted from submillimctcr/radio flux ratio analysis. We 
detected two lensed SMA sources toward A2390-3. The 
answer to the question of whether these two sources are 
physically related or the consequence of lensing distor- 
tion is still unclear. The fact that no counterparts are 
found in optical and infrared images indicates a heav- 
ily dust-embedded nature, and a lack of radio emission 
implies that both sources are located beyond z = 4.7. 



Our non-detection toward A2390-4 suggests either that 
there are multiple sources within the SCUBA beam or 
that the SCUBA detection may have been false. The 
accurate amplifications and fluxes from this work pro- 
vide a shallower slope in the faint-end 850 /im cumula- 
tive number counts. Our results emphasize the need for 
high-resolution observations toward single-dish detected 
SMGs. 

We gratefully acknowledge support from NSF grants 
AST 0709356 (C.C.C., L.L.C.) and AST 0708793 
(A.J.B.), National Science Council of Taiwan grant 99- 
2112-M-001-012-MY3 (W.-H.W.), the University of Wis- 
consin Research Committee with funds granted by the 
Wisconsin Alumni Research Foundation (A.J.B.), and 
the David and Lucile Packard Foundation (A.J.B.). 
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